Search for: All records

ABSTRACT High-velocity clouds (HVCs) may fuel future star formation in the Milky Way, but they must first survive their passage through the hot halo. While recent work has improved our understanding of the survival criterion for cloud–wind interactions, few observational comparisons exist that test this criterion. We therefore present an initial comparison of simulations with the Smith Cloud (SC; $d=$ 12.4 kpc, $$l, b = 40^{\circ }, -13^{\circ }$$) as mapped with the GALFA-HI (Galactic Arecibo L-Band Feed Array HI) survey. We use the SC’s observed properties to motivate simulations of comparable clouds in wind tunnel simulations with enzo-e, a magnetohydrodynamic code. For both observations and simulations, we generate moment maps, characterize turbulence through a projected first-order velocity structure function (VSF), and do the same for H i column density with a normalized autocovariance function. We explore how initial cloud conditions (such as radius, metallicity, thermal pressure, viewing angle, and distance) affect these statistics, demonstrating that the small-scale VSF is sensitive to cloud turbulence, while large scales depend on cloud bulk velocity and viewing angle. We find that some simulations reproduce key observational features (particularly the correlation between column density and velocity dispersion) but none match all observational probes at the same time (the large scales of the column density autocovariance is particularly challenging). We find that the simulated cloud (cloud C) showing growth via a turbulent radiative mixing layer (TRML) is the best match, implying the importance of TRML-mediated cooling for Milky Way HVCs. We conclude by suggesting improvements for simulations to better match observed HVCs. 

Abstract Turbulent radiative mixing layers play an important role in many astrophysical contexts where cool (≲10⁴K) clouds interact with hot flows (e.g., galactic winds, high-velocity clouds, infalling satellites in halos and clusters). The fate of these clouds (as well as many of their observable properties) is dictated by the competition between turbulence and radiative cooling; however, turbulence in these multiphase flows remains poorly understood. We have investigated the emergent turbulence arising in the interaction between clouds and supersonic winds in hydrodynamicenzo-esimulations. In order to obtain robust results, we employed multiple metrics to characterize the turbulent velocity,v_turb. We find four primary results when cooling is sufficient for cloud survival. First,v_turbmanifests clear temperature dependence. Initially,v_turbroughly matches the scaling of sound speed on temperature. In gas hotter than the temperature where cooling peaks, this dependence weakens with time untilv_turbis constant. Second, the relative velocity between the cloud and wind initially drives rapid growth ofv_turb. As it drops (from entrainment),v_turbstarts to decay before it stabilizes at roughly half its maximum. At late times, cooling flows appear to support turbulence. Third, the magnitude ofv_turbscales with the ratio between the hot phase sound-crossing time and the minimum cooling time. Finally, we find tentative evidence for a length scale associated with resolving turbulence. Underresolving this scale may cause violent shattering and affect the cloud’s large-scale morphological properties. 

Abstract We present a new suite of numerical simulations of the star-forming interstellar medium (ISM) in galactic disks using the TIGRESS-NCR framework. Distinctive aspects of our simulation suite are (1) sophisticated and comprehensive numerical treatments of essential physical processes including magnetohydrodynamics, self-gravity, and galactic differential rotation, as well as photochemistry, cooling, and heating coupled with direct ray-tracing UV radiation transfer and resolved supernova feedback and (2) wide parameter coverage including the variation in metallicity over $Z\prime \equiv Z/Z_{\odot }\sim 0.1-3$, gas surface density Σ_gas∼ 5–150M_⊙pc⁻², and stellar surface density Σ_star∼ 1–50M_⊙pc⁻². The range of emergent star formation rate surface density is Σ_SFR∼ 10⁻⁴–0.5M_⊙kpc⁻²yr⁻¹, and ISM total midplane pressure isP_tot/k_B= 10³–10⁶cm⁻³K, withP_totequal to the ISM weight $\mathcal{W}$. For given Σ_gasand Σ_star, we find ${\mathrm{\Sigma }}_{\mathrm{SFR}}\propto Z{\prime }^{0.3}$. We provide an interpretation based on the pressure-regulated feedback-modulated (PRFM) star formation theory. The total midplane pressure consists of thermal, turbulent, and magnetic stresses. We characterize feedback modulation in terms of the yield ϒ, defined as the ratio of each stress to Σ_SFR. The thermal feedback yield varies sensitively with both weight and metallicity as ${\mathrm{\Upsilon }}_{\mathrm{th}}\propto {\mathcal{W}}^{-0.46}Z{\prime }^{-0.53}$, while the combined turbulent and magnetic feedback yield shows weaker dependence ${\mathrm{\Upsilon }}_{\mathrm{turb}+\mathrm{mag}}\propto {\mathcal{W}}^{-0.22}Z{\prime }^{-0.18}$. The reduction in Σ_SFRat low metallicity is due mainly to enhanced thermal feedback yield, resulting from reduced attenuation of UV radiation. With the metallicity-dependent calibrations we provide, PRFM theory can be used for a new subgrid star formation prescription in cosmological simulations where the ISM is unresolved. 

ABSTRACT The origin of the cold phase in the circumgalactic medium (CGM) is a highly debated question. We investigate the contribution of satellite galaxies to the cold gas budget in the CGM of a Milky Way-like host galaxy. We perform controlled experiments with three different satellite mass distributions and identify several mechanisms by which satellites can add cold gas to the CGM, including ram pressure stripping and induced cooling in the mixing layer of the stripped cold gas. These two mechanisms contribute a comparable amount of cold gas to the host CGM. We find that the less massive satellites (≤109M⊙) not only lose all of their cold gas in a short period (∼ 0.5–1 Gyr), but their stripped cold clouds also mix with the hot CGM gas and get heated up quickly. However, stellar feedback from these less massive satellites can hugely alter the fate of their stripped gas. Feedback speeds up the destruction of the stripped cold clouds from these satellites by making them more diffuse with more surface area. On the other hand, the more massive satellites (LMC or SMC-like ∼1010M⊙) can add cold gas to the total gas budget of the host CGM for several Gyr. 

ABSTRACT In the absence of supplementary heat, the radiative cooling of halo gas around massive galaxies (Milky Way mass and above) leads to an excess of cold gas or stars beyond observed levels. Active galactic nucleus jet-induced heating is likely essential, but the specific properties of the jets remain unclear. Our previous work concludes from simulations of a halo with $$10^{14} \,\mathrm{ M}_\odot$$ that a successful jet model should have an energy flux comparable to the free-fall energy flux at the cooling radius and should inflate a sufficiently wide cocoon with a long enough cooling time. In this paper, we investigate three jet modes with constant fluxes satisfying the criteria, including high-temperature thermal jets, cosmic ray (CR)-dominant jets, and widely precessing kinetic jets in $$10^{12}-10^{15}\, {\rm M}_{\odot }$$ haloes using high-resolution, non-cosmological magnetohydrodynamic simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, conduction, and viscosity. We find that scaling the jet energy according to the free-fall energy at the cooling radius can successfully suppress the cooling flows and quench galaxies without violating observational constraints. On the contrary, if we scale the energy flux based on the total cooling rate within the cooling radius, strong interstellar medium cooling dominates this scaling, resulting in a jet flux exceeding what is needed. Among the three jet types, the CR-dominant jet is most effective in suppressing cooling flows across all surveyed halo masses due to enhanced CR pressure support. We confirm that the criteria for a successful jet model work across a wider range, encompassing halo masses of $$10^{12}-10^{15} {\rm M_\odot }$$. 

Abstract We make an in-depth analysis of different active galactic nuclei (AGN) jet models’ signatures, inducing quiescence in galaxies with a halo mass of 10¹²M_⊙. Three jet models, including cosmic-ray-dominant, hot thermal, and precessing kinetic jets, are studied at two energy flux levels each, compared to a jet-free, stellar feedback-only simulation. Each of our simulations is idealized isolated galaxy simulations with AGN jet powers that are constant in time and generated using GIZMO and with FIRE stellar feedback. We examine the distribution of Mgii, Ovi, and Oviiiions, alongside gas temperature and density profiles. Low-energy ions, like Mgii, concentrate in the interstellar medium (ISM), while higher energy ions, e.g., Oviii, prevail at the AGN jet cocoon’s edge. High-energy flux jets display an isotropic ion distribution with lower overall density. High-energy thermal or cosmic-ray jets pressurize at smaller radii, significantly suppressing core density. The cosmic-ray jet provides extra pressure support, extending cool and warm gas distribution. A break in the ion-to-mass ratio slope in Oviand Oviiiis demonstrated in the ISM-to-circumgalactic medium (CGM) transition (between 10 and 30 kpc), growing smoothly toward the CGM at greater distances. 

Abstract High-velocity outflows are ubiquitous in compact, massive (M_*∼ 10¹¹M_⊙),z∼ 0.5 galaxies with extreme star formation surface densities (Σ_SFR∼ 2000M_⊙yr⁻¹kpc⁻²). We have previously detected and characterized these outflows using Mgiiabsorption lines. To probe their full extent, we present Keck/KCWI integral field spectroscopy of the [Oii] and Mgiiemission nebulae surrounding all of the 12 galaxies in this study. We find that [Oii] is more effective than Mgiiin tracing low surface brightness, extended emission in these galaxies. The [Oii] nebulae are spatially extended beyond the stars, with radial extentR₉₀between 10 and 40 kpc. The nebulae exhibit nongravitational motions, indicating galactic outflows with maximum blueshifted velocities ranging from −335 to −1920 km s⁻¹. The outflow kinematics correlate with the bursty star formation histories of these galaxies. Galaxies with the most recent bursts of star formation (within the last <3 Myr) exhibit the highest central velocity dispersions (σ≳ 400 km s⁻¹), while the oldest bursts have the lowest-velocity outflows. Many galaxies exhibit both high-velocity cores and more extended, slower-moving gas indicative of multiple outflow episodes. The slower, larger outflows occurred earlier and have decelerated as they propagate into the circumgalactic medium and mix on timescales ≳50 Myr. 

Abstract The processes controlling the complex clump structure, phase distribution, and magnetic field geometry that develop across a broad range of scales in the turbulent interstellar medium (ISM) remain unclear. Using unprecedentedly high-resolution 3D magnetohydrodynamic simulations of thermally unstable turbulent systems, we show that large current sheets unstable to plasmoid-mediated reconnection form regularly throughout the volume. The plasmoids form in three distinct environments: (i) within cold clumps, (ii) at the asymmetric interface of the cold and warm phases, and (iii) within the warm, volume-filling phase. We then show that the complex magnetothermal phase structure is characterized by a predominantly highly magnetized cold phase, but that regions of high magnetic curvature, which are the sites of reconnection, span a broad range in temperature. Furthermore, we show that thermal instabilities change the scale-dependent anisotropy of the turbulent magnetic field, reducing the increase in eddy elongation at smaller scales. Finally, we show that most of the mass is contained in one contiguous cold structure surrounded by smaller clumps that follow a scale-free mass distribution. These clumps tend to be highly elongated and exhibit a size versus internal velocity relation consistent with supersonic turbulence and a relative clump distance–velocity scaling consistent with subsonic motion. We discuss the striking similarity of cold plasmoids to observed tiny-scale atomic and ionized structures and H i fibers and consider how the presence of plasmoids will modify the motion of charged particles, thereby impacting cosmic-ray transport and thermal conduction in the ISM and other similar systems. 

ABSTRACT Observed accretion rates onto the Milky Way and other local spirals fall short of that required to sustain star formation for cosmological timescales. A potential avenue for this unseen accretion is a rotating inflow in the volume-filling hot phase ($$\sim 10^6\, {\rm K}$$) of the circumgalactic medium (CGM), as suggested by some cosmological simulations. Using hydrodynamic simulations and a new analytic solution valid in the slow-rotation limit, we show that a hot inflow spins up as it approaches the galaxy, while remaining hot, subsonic, and quasi-spherical. Within the radius of angular momentum support ($$\sim 15\, {\rm kpc}$$ for the Milky Way) the hot flow flattens into a disc geometry and then cools from $$\sim 10^6$$ to $$\sim 10^4\, {\rm K}$$ at the disc–halo interface. Cooling affects all hot gas, rather than just a subset of individual gas clouds, implying that accretion via hot inflows does not rely on local thermal instability in contrast with ‘precipitation’ models for galaxy accretion. Prior to cooling and accretion the inflow completes ≈tcool/tff radians of rotation, where tcool/tff is the cooling time to free-fall time ratio in hot gas immediately outside the galaxy. The ratio tcool/tff may thus govern the development of turbulence and enhancement of magnetic fields in gas accreting onto low-redshift spirals. We show that if rotating hot inflows are common in Milky-Way-size disc galaxies, as predicted, then signatures of the expected hot gas rotation profile should be observable with X-ray telescopes and fast radio burst surveys. 

ABSTRACT In this paper, we study the filamentary substructure of 3.3 $$\mu$$m polycyclic aromatic hydrocarbon (PAH) emission from JWST/NIRCam observations in the base of the M 82 star-burst driven wind. We identify plume-like substructure within the PAH emission with widths of $$\sim$$50 pc. Several of those plumes extend to the edge of the field-of-view, and thus are at least 200–300 pc in length. In this region of the outflow, the vast majority ($$\sim$$70 per cent) of PAH emission is associated with the plumes. We show that those structures contain smaller scale ‘clouds’ with widths that are $$\sim$$5–15 pc, and they are morphologically similar to the results of ‘cloud-crushing’ simulations. We estimate the cloud-crushing time-scales of $$\sim$$0.5–3 Myr, depending on assumptions. We show this time-scale is consistent with a picture in which these observed PAH clouds survived break-out from the disc rather than being destroyed by the hot wind. The PAH emission in both the mid-plane and the outflow is shown to tightly correlate with that of Pa $$\alpha$$ emission (from Hubble Space Telescope data), at the scale of both plumes and clouds, though the ratio of PAH-to-Pa $$\alpha$$ increases at further distances from the mid-plane. Finally, we show that the outflow PAH emission reaches a local minimum in regions of the M 82 wind that are bright in X-ray emission. Our results are consistent cold gas in galactic outflows being launched via hierarchically structured plumes, and those small scale clouds are more likely to survive the wind environment when collected into the larger plume structure.